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'j^. ALEPH data on inclusive D*^ production in e^e~ annihilation. We take the charm- 

Qh! quark fragmentation function to be of the form proposed by Peterson et al. and 
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D*^ photoproduction in ep collisions at HERA are reasonably well reproduced. 
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1 Introduction 

Recently, the HI [l|] and ZEUS [^ collaborations at HERA presented data on the differen- 
tial cross section d?a /dy\ah dpT of inclusive D*"^ production in low-Q^ ep collisions, equiv- 
alent to photoproduction. Here yiah and px are the rapidity and transverse momentum 
of the produced D*^ mesons in the laboratory frame, respectively. These measurements 
extended up to pt = 12 GeV. Another source of information on D*^ production comes 
from e~^e~ —>■ D*^ + X at the Z-boson resonance |P,^. In this process, two mechanisms 
of D*"^ production contribute with similar rates. The D*^ mesons are either produced 
by Z ^ cc decay and subsequent c/c -^ D*^ fragmentation, or by Z ^ bb decay with 
subsequent fragmentation of h/h quarks into B mesons, which weakly decay into D*"^ 
mesons. The latter two-step process is usually treated as a one-step fragmentation pro- 
cess h/h -^ D*^. With the aid of very efficient bottom-tagging methods, these two main 
sources of inclusive D*^ production were disentangled with high purity, and separate cross 
sections for the two production mechanisms were presented. This allows one to determine 
separate fragmentation functions (FF's) for c/c -^ D*^ and h/h -^ D*^. Owing to the 
factorization theorem, these can then be used to make quantitative predictions for D*^ 
production in other reactions such as 7p -^ D*"^ + X, which is being measured at HERA. 
Such a program, which constitutes a test of the universality of the D*^ FF's, was recently 
carried out by two groups |]§,|§]. 

These works are based on the so-called massless-charm scheme [0. In this scheme, the 
charm-quark mass rUc is neglected, except in the initial conditions for the parton density 
functions (PDF's) and the FF's of the charm quarks. This should be a reasonable ap- 
proximation for center-of-mass (CM) energies a/s ^ mc in e^e~ annihilation or transverse 
momenta Pt ^ ?^c in 7P scattering. In this approach, charm is considered to be one of 
the active flavors inside the initial photons and protons, in the same way as the lighter -u, 
d, and s quarks. Then, the coUinear singularities which correspond to the as ln(/i^/m^) 
terms, where /x is an appropriate factorization scale (/i = y/s and /x = p-r, respectively), 
in a scheme where the charm quark is treated as a massive particle and only three ac- 
tive flavors are taken into account are absorbed into the charm-quark PDF's and FF's. 
Instead of absorbing the logarithmic terms, one can start with mc = and absorb the 
initial- and final-state coUinear singularities as usual in the modified minimal-subtraction 
(MS) factorization scheme, based on dimensional regularization. These two possibilities 
to achieve a massless-charm result, i.e., to absorb the coUinear singularities either in the 
form of logarithms or in the form of 1/e poles via MS factorization, differ by finite terms, 
as was shown in Ref. for the case of e"'"e~ -^ c/c + X. These terms can be considered 
as perturbative FF's at the low initial scale jIq of order ttIc, which are evolved to higher 
scales with the usual Altarelli-Parisi (AP) equations 0. Along these lines, cross sections 



for pp |jTO[, 7P [0, and 77 [|T2[ collisions were calculated in next-to-leading order (NLO), 
describing this way the transformation of a massless charm quark into a massive charm 
quark via a perturbative QCD cascade. Later, this approach was extended by including 
a nonperturbative FF which describes the transition from the massive charm quark and 
antiquark to the D*"^ mesons P,p!3|. In our earlier works on D*"^ production [|6|,p!4|, in- 



stead of working with the perturbative FF's, we incorporated their effect by an equivalent 
modification of the pure MS scheme with massless quarks, which led to a particular kind 
of massive subtraction scheme. This change of scheme was restricted to the final state. 

After fitting the nonperturbative fragmentation component to the e^e~ data taken by 
ARGUS ig and OPAL d in Ref. § and to those taken by ALEPH § and OPAL § in 



Ref. 10, the cross section of inclusive D*"^ photoproduction in ep collisions as measured by 
HI im and ZEUS could be reasonably well described. In particular, not only the shape, 
but also the normalization of the cross section was very well accounted for 0. These 
papers showed that, in order to successfully predict inclusive D*"^ production in low-Q^ 
ep collisions, it is essential to incorporate information on the fragmentation process from 
other reactions and to use the very same factorization scheme for all considered processes. 

The fragmentation of charm or bottom quarks into D*^ mesons cannot be calculated 
from first principles in perturbative QCD. In fact, in order to realistically describe the 
formation of the D*^ mesons, a nonperturbative component, which is not known theo- 
retically, is always needed. Hence, it is certainly appropriate to give up the perturbative 
component of the FF input altogether and to describe the c/c -^ D*^ and b/b -^ D*^ 
transitions entirely by nonperturbative FF's, as is usually done for the fragmentation of 
M, d, and s quarks into light mesons. 

The aim of the present work is to reconsider the production of D*^ mesons both in e+e^ 
annihilation and photoproduction adopting this puristic approach of fully nonperturbative 
fragmentation in the MS scheme with five massless flavors. This analysis is based on very 
recent high-statistics data from the OPAL [T^ and ALEPH ||l^ collaborations at LEPl. 
The FF's thus fitted are then used to update our predictions for the cross section of 
D*"^ photoproduction in ep collisions to be compared with new data from the ZEUS 
collaboration HTB . 



The outline of this paper is as follows. In Sec. 2, we briefly recall the theoretical 
framework underlying the extraction of FF's from e^e~ data, which has already been 
introduced in Refs. [§,0. Then, we present the FF's we obtained by fitting the new D*^ 



data from OPAL |T6l and ALEPH ^ at leading order (LO) and NLO in the pure MS 
factorization scheme with five massless flavors. In Sec. 3, we apply the nonperturbative 
parameters determined in Sec. 2 to make LO and NLO predictions for D*^ photopro- 
duction in ep collisions, which we then compare with the latest ZEUS data |T^. Our 
conclusions are summarized in Sec. 4. 

2 D*^ Production in e^e~ Annihilation 

Our procedure to construct LO and NLO sets of D*^ FF's has already been described 
in our previous paper [01. Here, we only give those details which differ from Ref. H. As 



experimental input, we now use the new LEPl data from OPAL [jT6l and ALEPH |T] 
While the older ALEPH [^ and OPAL data |^ agreed well enough to be simultaneously 
fitted 0, the new data samples do not sufficiently overlap. Thus, we refrain from per- 
forming a combined fit. Instead, we generate independent LO and NLO FF sets for the 



two data samples. This allows us to check actually how well our predictions for other 
/)*±_pi-oduction cross sections are constrained by the available information from e~^e~ an- 
nihilation, an error source which has so far been rather difficult to assess quantitatively. 
In e"'"e~ annihilation at the Z-boson resonance, charmed mesons are produced either 
directly through the hadronization of charm quarks produced hj Z —>■ cc or via the 
weak decays of B hadrons from Z —>■ bb, with an approximately equal rate. Charmed 
mesons from Z ^ cc allow us to determine the charm-quark FF. The main achievement 
of both experimental analyses [16,171 was to disentangle these two main sources of D*^ 



production at the Z-boson resonance. In Ref. |T^, a combination of several charm-quark 



tagging methods, based on fully and partially reconstructed D*^ mesons, and a bottom- 



quark tag, based on identified electrons and muons, was used. In Ref. IT^], the separation 
of cc and bb events was improved by means of an algorithm based on the measurement 
of the impact parameter of charged tracks and the effective mass of those least fitted to 
the primary vertex. As in the earlier analyses 1^0, the D*^ mesons were identified via 
the decay chain D*~^ -^ D^tt^ and D° — > K^tc^ (and the analogous chain for the D*^ 
meson), which allows for a particularly clean signal reconstruction. The experimental 
cross sections |[T6|JT7|1 were presented as distributions differential in a; = 2E{D*'^)/y/s, 
where E{D*^) is the measured energy of the D*^ candidate, which are normalized to the 
total number of multihadronic Z-boson decays. Besides the total D*^ yield, which receives 
contributions from Z ^ cc and Z ^ bb decays as well as from gluon fragmentation, both 
experimental groups separately specified results for D*^ mesons from tagged Z —>■ bb 
events. The contribution due to charm-quark fragmentation is peaked at large x, whereas 
the one due to bottom-quark fragmentation has its maximum at small x. 

For the fits, we use the x bins in the interval [0.1, 1.0] and integrate the theoretical 
functions over the bin widths, which is equivalent to the experimental binning procedure. 
As in the experimental analyses, we sum over D*'^ and D*~ mesons. As a consequence, 
there is no difference between the FF's of a given quark and its antiquark. As in Ref. [P, 
we take the starting scales for the D*"^ FF's of the gluon and the u, d, s, and c quarks 
and antiquarks to be fiQ = 2mc, while we take /xq = 2rrib for the FF's of the b quark and 
antiquark. The FF's of the gluon and the first three fiavors are assumed to be zero at the 
starting scale. These FF's are generated through the /x^ evolution, and the FF's of the 
first three quarks and antiquarks coincide with each other at all scales /i. 

We employ two different forms for the parameterization of the charm- and bottom- 
quark FF's at their respective starting scales. In the case of charm, we use the Peterson 
distribution [jl9[. 

In the case of bottom, we adopt the ansatz 

Dt{x,fil) = Nx''{l-xf, (2) 

which is frequently used for the FF's of light hadrons. However, in contrast to the light- 
hadron case, we expect that a > for D*^ mesons, so that Db{x,fiQ) — > as a; -^ 0. 



Furthermore, we expect that a < /3, since Db{x,fil) is supposed to have its maximum 
somewhere at x < 0.5. This choice of starting distributions was also used for the mixed 
set (M) in Ref. p. The Peterson form is particularly suitable for FF's that peak at large 
X. Since the bottom-quark FF, being a convolution of the b ^ B fragmentation and the 
subsequent B —>■ D*^ + X decay, has its maximum at small x values, we obtain intolerably 
bad fits if we also use Eq. (|lD in this case. In Ref. 0, we also studied the alternative where 
Eq. (H) is used both for the charm and bottom quarks, leading to the standard set (S). 
Although set S yielded slightly lower x^ values per degree of freedom (xdf) ^^ ^^^ expense 
of having one more fit parameter, both sets, S and M, led to almost identical predictions 
for D*^ photoproduction in the kinematical regions of present interest at HERA. In this 
work, we thus limit ourselves to the mixed ansatz. This facilitates the comparison of 
our results, especially on ec, with the literature on charm-quark fragmentation, where the 
Peterson ansatz (|l]) is commonly used. Ansatz (^ was also employed some time ago in 
Ref. O] to describe the nonperturbative FF's of charm quarks into D and D* mesons 
and those of bottom quarks into B mesons. 

The calculation of the cross section {l/crtot)da/dx for e^e^ ^ •y, Z ^ D*^ + X is 
performed as described in Ref. 0, except that we now abandon the subtraction terms 
dQa{x) specified in Eq. (7) therein, for reasons explained in Sec. 1. All relevant formulas 
and references may be found in Ref. 0. 

Both OPAL [|l^ and ALEPH fl^ presented momentum distributions for their full 



D*^ samples and for their Z ^ bb subsamples. We received these data in numerical 
form via private communications [T^JT^. The OPAL data are displayed in Fig. 4 of 
Ref. [|I^ in the form {l/N^g^d)dNij*±/dx, where N£)*± is the number of D*^ candidates 
reconstructed through the decay chain mentioned above. In order to convert this into 
the cross section {l/atot)dcr/dx, we need to divide by the branching fractions B{D*^ -^ 
D^TT+) = 0.683 ± 0.014 and B{D° -^ R-n^) = 0.0383 ± 0.0012 ^. The momentum 
distribution of the full ALEPH sample, which is shown in Fig. 3 of Ref. fl^, has to be 
treated in the same way. In the case of the ALEPH Z -^bb subsample, which is presented 
in Fig. 8 of Ref. [0, we need to include the additional factor i?f, x /f, = 0.0464, where 
Rb = T[Z —>■ bb)/r[Z — >■ hadrons) and fb is the fraction of bb events in the sample which 
remains after subtracting the events from gluon splitting; this factor may be extracted 
from Ref. [|17|. The ALEPH data |T^ have slightly smaller statistical errors than the 



OPAL data 16 



As for the asymptotic scale parameter appropriate for five active quark flavors, we 
adopt the LO (NLO) value A^ = 108 MeV (227 MeV) from our study of inclusive 



charged-pion and -kaon production |21|. The particular choice of A^ is not essential; 



MS 
(5) 

other values, in particular higher LO values of A|-t^, can easily be accommodated by 
slightly shifting the other fit parameters without changing the quality of the fit. As 
in Ref. P], we take the charm- and bottom-quark masses to be rric = 1.5 GeV and 
nib = 5 GeV, respectively. Since rric and m& only enter via the definitions of the starting 
scales fiQ of the FF's, their precise values are immaterial for our fit. 

The values of N and ec in Eq. (|I]) and of A^, a, and /5 in Eq. (||) which result from 



our LO and NLO fits to the OPAL and ALEPH data are summarized in Table |1|. In the 
following, we refer to the corresponding FF's as sets LO O, NLO O, LO A, and NLO A, 
respectively. In Table 0, we list three Xdf values for each of these four fits: one for the 
Z — > 66 subsample; one for the total sample (sum of cc-tagged, 66-tagged, and gluon- 
splitting events); and the average evaluated by taking into account both the Z -^ bb 
subsample and the total sample. We observe that all four fits are quite successful, with 
Xqp values of order unity. The LO and NLO fits to the OPAL data are slightly better 
than those to the ALEPH data. In each case, the Z — > 66 subsample tends to be less well 
described by the fit than the total sample. As expected on general grounds, the NLO fits 
are superior to the LO fits. 

In Fig. |l]a, we compare the OPAL data |jl6|] with the LO and NLO calculations using 



sets LO O and NLO O, respectively. The analogous analysis for the ALEPH data |]T^ is 
shown in Fig. [l|b. Notice that the distributions plotted in Figs. |l]a and b correspond to 
{l/atot)da/dx, i.e., the experimental data are multiplied by the factors specified above. 
Except at very small x, the LO and NLO results are very similar. This is also true for 
the distributions at the starting scale, as may be seen by comparing the corresponding 
LO and NLO parameters in Table ^ Only ec changes appreciably as we pass from LO 
to NLO. The branching of the LO and NLO results at small x indicates that, in this 
region, the perturbative treatment ceases to be valid. This is related to the phase-space 
boundary for the production of D*^ mesons at Xmin = '2'fn{D*^)/-s/s, where m{D*^) is the 
D*"^ mass. At ^/s = Mz, one has Xmin = 0.046. This is approximately where our NLO 
results turn negative. Since our massless- quark approach is not expected to be valid in 
regions of phase space where f\.mte-m{D*^) effects are important, our results should only 
be considered meaningful for x^Xcut = 0.1, say. We also encountered a similar small-x 
behaviour in our previous analysis [^, where we fitted the older ALEPH [^] and OPAL 
data in the framework of the massive subtraction scheme. Since the two rightmost 
data points of the ALEPH Z —>■ bb sample, corresponding to the x bins [0.85,0.9] and 
[0.9,0.95], come with negative cross sections and rather small errors |]T^, we excluded 
them from our fits. 

Comparing Figs. ^ and b, we observe that the OPAL and ALEPH data are indeed 
somewhat different in shape and do not mutually overlap within their errors. Thus, 
combined LO and NLO fits would lead to intolerable values of Xdf- By the same token, 
the separate fits to the OPAL and ALEPH data lead to significantly different values for 
the parameters appearing in the starting distributions (|l|) and (|]), as may be seen from 
Table 0. There are also striking differences between the parameters resulting from the LO 
and NLO fits to the same data sets. The value ec = 0.0851 of set LO O is very similar 
to the value ec = 0.0856 of set LO M in our previous analysis [^]. On the other hand, 
our new NLO results for ec, namely 0.116 for set NLO O and 0.185 for set NLO A, are 
significantly larger than our previous value 0.0204 for set NLO M |^. This dramatic shift 
in ec is due to the fact that we are now using a different scheme for the factorization 
of the final-state collinear singularities (namely, the pure MS scheme with five massless 
flavors) than in Ref. ^, where we modified this scheme so as to incorporate the effect of 
the perturbative FF's ||^ (massive subtraction scheme). One important advantage of our 
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new approach becomes apparent if we compare Figs. |I|a and b with Figs, lb and 2b in 
Ref. [1; see also Fig. 1 in Ref. §. While in the NLO analyses of Refs. |J|, the NLO 
cross section of e^e~ — *■ D*"^ + X turned negative for x^O.9, the new NLO calculation 
stays positive in the upper x range, coincides there with the LO calculation, and nicely 
describes the data. We attribute the unphysical large-x behaviour of the NLO calculation 
in Refs. 0J^ to the use of the massive subtraction scheme in conjunction with the special 
forms for the nonperturbative FF's at the starting scale. This led to the appearance of 
large Sudakov logarithms which spoiled the NLO result at a; ^0.9. This feature forced us 
P], and also the authors of Ref. p, to exclude from the fits the rather precise experimental 
data points at x close to unity. By contrast, the NLO calculation in the pure MS scheme 
does not suffer from theoretical problems at high x, and there is no need to omit valuable 
experimental information in that region. 

As mentioned above, we take the FF's of the partons g, u, it, d, d, s, s to be vanishing 
at their starting scale /iq = 2mc. However, these FF's are generated via the AP evolution 
to the high scale /x = ^/s. Thus, apart from the FF's of the heavy quarks c, c, b, b, also 
these radiatively generated FF's contribute to the cross section. All these contributions 
are properly included in the total result for {l/atot)do'/dx shown in Figs. |l|a and b. At 
LEPl energies, the contribution from the first three quark flavors is still negligible; it is 
concentrated at small x, and only amounts to a few percent of the integrated cross section. 
However, the contribution from the gluon FF, which appears at NLO in connection with 
three-parton final states, is numerically significant. Motivated by the decomposition of 
{l/ax,ot)da/dx in terms of parton-level cross sections (see Eq. (5) in Ref. 0), in Figs. |l]a 
and b, we distributed this contribution over the Z ^ cc and Z ^ bb channels in the 
proportion e^ : e^, where Cg is the effective electroweak coupling of the quark q to the Z 
boson and the photon including propagator adjustments. This procedure approximately 
produces the quantities that should be compared with the OPAL and ALEPH data. We 
ignore the effect of electromagnetic initial-state radiation, which has not been corrected 
for in the data. 

As in Ref. Q, we study the branching fractions for the transitions of charm and bottom 
quarks to D*"^ mesons, defined by 



Bgifi) = / dxDQ{x,fi ), (3) 

where Q = c,b and Xcut = 0.1. This allows us to test the consistency of our fits with 
information presented in the experimental papers ||T6|JT7[| . The contribution from the 
omitted region < a; < Xcut is close to zero. For all four FF sets, we calculate Bq{^) 
at the respective threshold /x = 2mQ and at the Z-boson resonance /x = Mz and present 
the outcome in Table |^. As expected, the values of Bqi^fi) change very little under the 
evolution from 2mQ to Mz, and they are very similar for Q = c,b. The OPAL analysis ||T6| , 



which is conceptually very different from ours, yielded Bc{Mz) = 0.222±0.014±0.014 and 
Bh{Mz) = 0.173 ± 0.016 ± 0.012, where the first (second) error is statistical (systematic). 
These values lie in the same ball park as our corresponding LO-0 and NLO-0 results 
in Table ^ The ALEPH paper |1^ does not explicitly quote values for these branching 



fractions. However, by combining other results given in Section 8 of Ref. [^, we infer 
that Bc{Mz) = 0.230 ^QQgg, which nicely agrees with our corresponding LO-A and NLO-A 
results in Table ^. 

Another quantity of interest, which can directly be compared with experiment, is the 
mean momentum fraction, 

1 /•! 
{x)q{^i) = -——- dxxDQ{x,fi), (4) 

where Q = c,b. In Table H, we collect the values of {x)Q{fi) for Q = c,b evaluated at 
fi = 2mQ, Mz with the four FF sets. At fixed fi, the differences between the OPAL and 
ALEPH sets and between LO and NLO are marginal. However, the AP evolution from 
fi = 2mQ to yU = Mz leads to a significant reduction of {x)q{^), especially in the case 
oi Q = c. Our results for {x)c{Mz) should be compared with the experimental numbers 
reported by OPAL [0] and ALEPH 0, 



{x)c{Mz) = 0.515 ±0.002 ±0.009 (OPAL), 

{x)c{Mz) = 0.4878 ±0.0046 ±0.0061 (ALEPH), (5) 

respectively. The central values are very close to those reported in the previous OPAL 
0] and ALEPH publications, but the errors are now considerably smaller. Comparing 
Eq. (|^) with Table |, we conclude that our results for {x)c{Mz) slightly undershoot the 
independent experimental determinations in Refs. |lT6| , p!7| , by about 5% (11%) at LO 



(NLO). In order to assess the theoretical uncertainty of the results in Table ^ related to 
the choice of ansatz at the starting scale, we recall that, in Ref. 0, the LO-S (NLO-S) 
result for {x)c{Mz) turned out to be 6% (4%) larger than the LO-M (NLO-M) result. At 
this point, we should also mention that the results in Eq. (^ are not directly extracted 
from the measured x distributions, but from calculations based on some Monte-Carlo 
model with parameters fitted to the experimental data. Whether this is the actual source 
of the difference remains unclear for the time being. 

Having constructed D*^ FF's from LEPl data, it is interesting to quantitatively inves- 
tigate whether they lead to a consistent description of available data on e^e~ -^ D*^ + X 
at other values of ^/s. This would represent a direct test of the underlying scaling viola- 
tion of fragmentation as implemented in the QCD-improved parton model via the timelike 
AP equations. To this end, we select the data from ARGUS [|T^] at ^/s = 10.49 GeV, 
from HRS |2| at v^ = 29 GeV, and from TASSO p3| at v^ = 34.2 GeV. In Fig. |a, we 



compare these data with our respective LO and NLO predictions based on sets LO O and 
NLO O, respectively. For reference, also the OPAL data |T^ are included. The analogous 
comparison based on sets LO A and NLO A is shown in Fig. |]b. To quantitatively assess 
these comparisons, we summarize the corresponding values of Xdf i^ Table ^. We observe 
that the scaling violation encoded in the experimental data is faithfully described by our 
theoretical predictions. As expected on general grounds, the Xdf values tend to be lower 
for the NLO analyses. Within the framework of the QCD-improved parton model, the 
OPAL data appear to be more consistent with the ARGUS data than the ALEPH data, 
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while the ALEPH data seem agree better with the HRS data than the OPAL data do. 
However, we should bear in mind that, in contrast to the LEPl data, the ARGUS, HRS, 
and TASSO data are presented in the form sda/dx and thus suffer from an additional 
normalization uncertainty. 

The successful comparisons in Figs. 0a and b and Table |^ reassure us that our D*^ 
FF's, although constructed at ^/s = Mz, also lead to useful descriptions of D*^ frag- 
mentation at other scales. In the next section, we exploit this property together with the 
universality of fragmentation to make predictions for inclusive D*^ photoproduction at 
HERA. 

3 D*^ Production in Low-Q^ ep Collisions 

In this section, we compare our NLO predictions for the cross section of inclusive D*^ pho- 
toproduction in ep scattering at HERA with the 1996 data from the ZEUS collaboration, 
which were presented by Y. Eisenberg at the 1997 International Europhysics Conference 
on High Energy Physics [|18|. We emphasize that these data are still preliminary. 

The present HERA conditions are such that Ep = 820 GeV protons collide with 
Ee = 27.5 GeV positrons in the laboratory frame. The rapidity is taken to be positive in 
the proton flight direction. The quasi-real photon spectrum is described in the Weizsacker- 
Williams approximation by Eq. (5) of Ref. [0. This spectrum depends on the photon- 
energy fraction, x = E^/E^, and the maximum photon-virtuality, Q^ax- ^^ the ZEUS 
experiment @,0, where the final-state electron is not detected, Q^ax — ^ GeV^ and 
0.147 < X < 0.869, which corresponds to 7p CM energies in the range 115 GeV < 
W < 280 GeV. We adopt all these kinematic conditions in our analysis. We work at 
NLO in the MS scheme with Uf = 4 flavors. For the proton and photon PDF's we 
use set CTEQ4M [|3 with A^^ = 296 MeV and set GRV HO |2| converted to the 
MS factorization scheme, respectively. We evaluate ^^(/i^) from the two-loop formula 
with this value of A^^. The A^J- values implemented in the photon PDF's and the 

D*^ FF's are 200 MeV and 352 MeV, which corresponds to the value A^ = 227 MeV 
quoted in Sec. 2, respectively. We identify the factorization scales associated with the 
proton, photon, and D*^ mesons, and collectively denote them by Mf. We choose the 
renormalization and factorization scales to be /i = rrij' and Mf = 2mT, respectively. 



where rriT = yPr + ^ is the D*"^ transverse mass. Whenever we present LO results, 
these are consistently computed using set CTEQ4L [Bl] of proton PDF's, set GRV LO 



| 25| of photon PDF's, the LO versions of our two sets of D* FF's, the one-loop formula 



for as with AjJ- = 236 MeV |2^ , and the LO hard-scattering cross sections. 
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The photoproduction cross section is a superposition of the direct- and resolved-photon 
contributions. In our NLO analysis, the resolved-photon contribution is larger than the 
direct one for moderate px- This statement depends, however, on the factorization scheme; 
only the sum of both contributions is a physical observable and can be compared with 
experimental data. The bulk of the resolved-photon cross section is due to the charm 



content of the photon []I4 . 



The ZEUS data [|T8l come as three distributions: (i) the px distribution da/dpT in- 
tegrated over —1.5 < yiah < 1 and 115 GeV < W < 280 GeV; (ii) the y\ah distribution 
da/dyub integrated over 4 GeV < pt < 12 GeV and 115 GeV < W < 280 GeV; and (iii) 
the W distribution da/dW integrated over 4 GeV < pr < 12 GeV and — 1.5 < yiab < 1- 
In Figs. I^a-c, we compare the measured distributions (i)-(iii) with our LO and NLO pre- 
dictions based on FF sets LO O, LO A, NLO O, and NLO A. In all three cases, the LO-0 
(NLO-0) result slightly exceeds the LO-A (NLO-A) result, but the differences are still 
small compared to the experimental uncertainty of the ZEUS data. This means that the 
details of charm fragmentation are sufficiently well constrained by the LEPl data, and 
that the HERA data are not yet precise enough to resolve the presently existing differ- 
ence between the OPAL [^ and ALEPH |T^ data. In Figs, ^-c, the shifts due to use of 



different experimental input (OPAL versus ALEPH) tend to be less significant than the 
effects due to the inclusion of higher-order corrections (LO versus NLO). However, the 
differences between the LO and NLO results are still considerably smaller than the errors 
on the ZEUS data, and they are approximately equal for the OPAL and ALEPH sets. The 
similarity of the LO and NLO results indicates good perturbative stability. Apparently, 
this is a special virtue of the pure MS factorization scheme with strictly massless quark 
flavors, which we are using here. By contrast, we encountered a sizeable gap between the 
LO and NLO results in the massive subtraction scheme 0. 

In Fig. I^a, the NLO distributions fall off slightly less strongly with increasing px than 
the LO ones, and thus agree better with the data at large px- The overall agreement 
with the data is very satisfactory, even at small pt, where the massless approach ceases 
to be valid. The theoretical predictions somewhat undershoot the measurement in the 
highest Pt bin, but the experimental error is still rather sizeable there. As for Fig. ^, 
the effect of the higher-order corrections and the influence of the experimental input is 
most pronounced in the backward direction, where all theoretical predictions agree with 
the experimental data within one standard deviation. At yiab = —1.5, the NLO-0 result 
exceeds the LO-0 one by 25% and the NLO-A one by 15%. On the other hand, in the 
forward direction, at j/iab > 0, the ZEUS data points signiflcantly overshoot the theoretical 
predictions. There, the experimental errors are largest, while the theoretical uncertainty 
is relatively small. In fact, both the NLO corrections and the uncertainty due to spread 
of the LEPl data are negligibly small there. Furthermore, we know from Ref. [^ that the 
speciflc choice of ansatz for the charm FF's at the starting scale is numerically irrelevant. 
Changes in Aj^ or ^q are expected to be compensated through the flt by appropriate 
shifts in the input parameters of Eqs. (|l]) and (||), and do not represent an appreciable 
theoretical error source either. The variation in cross section due to the use of different 



proton PDF's is just of order 10%, as long as up-to-date sets are considered [0]. The most 



substantial source of theoretical uncertainty is related to the photon PDF's. Among the 



different sets used in Ref. |]§], GRV HO |^ gave the smallest cross section at yiab > 0; the 
variation was found to be about 40%. An increase of this size would render the theoretical 
predictions compatible with the ZEUS data point at yiab = 0.75, while it would be far to 
small to explain the data point at yiab = 1-25, which represents an unexpected rise in cross 
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section. It remains to be seen if such a rise will also be observed by the HI collaboration. 
Looking at Fig. ^c, we observe that, within the rather large experimental errors, all four 
theoretical predictions agree reasonably well with the data. The NLO-0 result leads to 
the best description of the data. Only the data point at H^ = 252.5 GeV is significantly 
above the theoretical expectation. 



In the ZEUS paper [|18|, the j/iab and W distributions were also presented for a 
minimum-pr cut of 3 GeV. We performed similar comparisons with these data, too, but 
do not display them here, since our massless- charm and -bottom scheme is more appro- 
priate for high values of pr- In the case of the yiab distribution, the general picture is 
very similar to Fig. ^d, except that the data point at yiab = 0.75 is then consistent with 
the theoretical expectation. Similarly, in the case of the W distribution, the data point 
at W = 252.5 GeV then agrees considerably better with the theoretical prediction. 

4 Conclusions 

The OPAL and ALEPH collaborations presented measurements of the fractional energy 
spectrum of inclusive D*^ production in Z-boson decays based on their entire LEPl data 
samples |]16| , p^ . Apart from the full cross section of D*^ production, they also determined 
the contribution arising from Z -^ bb decays. This enabled us to update our LO and NLO 
fits of the D*^ FF's 0, which were based on previous OPAL and ALEPH analyses with 
considerably lower statistics 0,01. 

At the same time, we also incorporated a conceptual modification of our theoretical 
framework, so as to eliminate two minor weaknesses of our previous approach 0: (i) 
At NLO, the differential cross section da/dx of e^e" -^ D*^ + X turned negative for 
x^O.9, so that the experimental data points in this region needed to be excluded from 
the fits; this problem was also encountered in Ref. [^. (ii) The LO and NLO results for 
the cross section of inclusive D*^ photoproduction in ep collisions substantially differed 
from each other, indicating that, in this particular scheme, higher orders beyond present 
control are likely to be significant. In Ref. 0, the factorization of the final-state coUinear 
singularities associated with the charm and bottom quarks was performed in a scheme 
which differs from the pure MS scheme with five massless quark flavors by the subtrac- 
tion of certain finite functions from the parton-level cross sections. These functions were 
chosen in such a way that the resulting change of scheme is equivalent to the use of the 
perturbative FF's of Ref. [§. Since these perturbative FF's describe the production of 
a massive quark Q rather than a heavy hadron if, they need to be complemented by a 
nonperturbative component which accounts for the Q —>■ H fragmentation process. In 
particular, this is requisite to make sure that the Q ^ H branching ratio and the average 
H to Q longitudinal-momentum fraction take their measured values, which are indeed 
smaller than unity. Since the charm quark is only moderately heavy, the effect of this 
nonperturbative component is rather dramatic, and its omission renders the theoretical 
description completely inadequate, as was demonstrated in Ref. [|14|. Clearly, the convo- 



lution of a perturbative FF with a nonperturbative ansatz can in turn be considered as a 



11 



nonperturbative FF. This consideration brings us back to the pure MS scheme with five 
massless flavors, from which we started off. In this paper, we advocated the apphcation 
of this scheme to charmed-meson production at high ^/s in e~^e~ annihilation and at large 
Pt in ep scattering or similar types of collision. As a matter of principle, the question 
whether this approach is justified in the present case or not can only be answered by 
nature itself. We found that the pure MS scheme is clearly favoured by a wealth of exper- 
imental data on D*^ production |[T5|JT6|JT7|JT8| , |22| , ^ . In fact, both drawbacks mentioned 



above are nicely avoided. Furthermore, comparing Figs. |^a and b with Fig. 3 in Ref . ^ , 
we observe that the general description of the e+e~ data from DORIS [|l^, PEP p2|, and 
PETRA p3l is improved, especially in the upper x range. 



It is important to bear in mind that the fit results for the input parameters in Eqs. (|I|) 
and (0), including the value of Peterson's ec, are highly scheme dependent at NLO, and 
must not be naively compared disregarding the theoretical framework which they refer to. 
If we compare the values of N and ec for the charm FF of set NLO M previously obtained 
in the massive subtraction scheme with those of set NLO O appropriate to the pure 
MS scheme, we find dramatic differences: N = 0.0677 and e^ = 0.0204 for set NLO M 
versus A^ = 0.267 and ec = 0.116 for set NLO O. On the other hand, the respective LO 
results are almost identical (A^ = 0.202 and e^ = 0.0856 for set LO M versus A^ = 0.223 
and ec = 0.0851 for set LO O), which indicates that the underlying data (from Refs. 0J^ 
for set M and from Ref. ||16[ for set O) are consistent with each other. 

Finally, we should caution the reader that the massless-quark approximation used 
here and the resulting FF's are only appropriate for processes of D*^ production which 
are characterized by an energy scale that is large against the charm-quark mass. In 
particular, this approach should not be expected to yield meaningful predictions for the 
photoproduction of small-py D*^ mesons at HERA. Therefore, in order to substantiate 
the comparison of LEPl and HERA data in the framework of the QCD-improved parton 
model endowed with nonperturbative FF's, it would be very desirable if the statistics of 
the HERA data was increased in the large-pr bins. 
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TABLE CAPTIONS 

Table 1: Fit parameters for the charm- and bottom-quark FF's of sets LO O, LO A, 
NLO O, and NLO A. The corresponding starting scales are /io = ^rric = 3 GeV and 
/io = 2mh = 10 GeV, respectively. All other FF's are taken to be zero at fio = 27?2c 



Table 2: x P^i' degree of freedom pertaining to the LO and NLO fits to the OPAL [16 



and ALEPH \T7\ data. In each case, Xdf is calculated for the Z —>■ bb sample, the full 
sample, and the combination of both. 

Table 3: Branching fractions of charm and bottom quarks into D*^ mesons at the re- 
spective starting scales and at /i = Mz evaluated from Eq. (^ with sets LO O, LO A, 
NLO O, and NLO A. 

Table 4: Average momentum fractions of D*^ mesons produced through charm- and 
bottom-quark fragmentation at the respective starting scales and at /i = Mz evaluated 
from Eq. (|) with sets LO O, LO A, NLO O, and NLO A. 

Table 5: x^ P^i' degree of freedom evaluated with sets LO O, LO A, NLO O, and NLO A 
for the ARGUS fg, HRS ^, and TASSO |23[ data. TASSOl (TASSO 2) refers to the 



D° -^ K-n+n+n- {D° -^ K-TX+) channel. 

FIGURE CAPTIONS 

Figure 1: (a) The cross sections of inclusive D*"^ production in e^e~ annihilation evaluated 
with sets LO O and NLO O are compared with the OPAL data [0. The three pairs of 
curves correspond to the Z — > cc, Z — >^ 66, and full samples, (b) The same for sets LO A 



and NLO A and the ALEPH data WJ 



Figure 2: (a) The ARGUS H, HRS H, TASSO [H, and OPAL pf data on inclusive 
D*"^ production in e^e~ annihilation are compared with the predictions based on sets 



LO O and NLO O. (b) The same for sets LO A and NLO A and the ALEPH data 1T7| 



For separation, the data have been rescaled by powers of 10. In the case of TASSO, the 
open triangles refer to the D^ — >■ K~tt^'k^tt~ channel (TASSO 1) and the solid triangles 
to the D^ -^ K-TI+ channel (TASSO 2). In the cases of OPAL and ALEPH, we consider 
the dimensionless quantity {l/atot)dcr/dx. 

Figure 3: The predictions of inclusive D*^ photoproduction in ep collisions based on sets 
LO O, LO A, NLO O, and NLO A are compared with the ZEUS data [|l|]. We consider (a) 
thepT distribution da/dpT integrated over —1.5 < yiab < 1 and 115 GeV < W < 280 GeV, 
(b) the yiab distribution da/dyiah integrated over 4 GeV < pr < 12 GeV and 115 GeV < 
W < 280 GeV, and (c) the W distribution da/dW integrated over 4 GeV < pr < 12 GeV 
and -1.5 < yi^h < 1- 
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